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Abstract 


Several fiber bridging models were reviewed and applied in this research to study the matrix 
fatigue crack growth behavior in center notched [0]^ SCS-6/TM5-3 and [0]^ SCS-6/Ti-6Al-4V 

laminates. Observations revealed that fatigue damage consisted primarily of matrix cracks and fiber- 
matrix interfacial failure in the [0]g SCS-6/TM5-3 laminates. Fiber-matrix interface failure included 

fracture of the brittle reaction zone and cracking between the two carbon rich fiber coatings. Intact fibers 

in the wake of the matrix cracks reduce the stress intensity factor range. Thus, an applied stress intensity 
factor range (AK a __ = AS^/rta ) is inappropriate to characterize matrix crack growth behavior. Fiber 

bridging models were used to determine the matrix stress intensity factor range in titanium metal matrix 
composites. In these models, the fibers in the wake of the crack are idealized as a closure pressure. An 
unknown constant frictional shear stress is assumed to act along the debond or slip length of the bridging 
fibers. In this study, the frictional shear stress was used as a curve fitting parameter to available data 
(crack growth data, crack opening displacement data, and debond length data). Large variations in the 
frictional shear stress required to fit the experimental data indicate that the fiber bridging models in their 
present form lack predictive capabilities. However, these models provide an efficient and relatively 
simple engineering method for conducting parametric studies of the matrix crack growth behavior based 
on constituent properties. 
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Nomenclature 

Current crack length and unbridged initial crack length, mm 
Superscript referring to crack centerline 

2 

Cross sectional area of fiber and matrix, m 

Superscript referring to end line of slip region 

Paris crack growth coefficient constants, yrn/MPa*cyclc 
Fiber and matrix modulus, MPa 

Composite longitudinal modulus, MPa 

Composite transverse modulus, MPa 

Composite shear modulus, MPa 

Mode I stress intensity factor range, MPaVm 
Applied stress intensity factor range, MPaVm 

Discrete stress intensity factor range in matrix, MPaVm 
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Continuum stress intensity factor range in composite, MPaVm 
Slip length, m 

Slip length of frictional shear stress reversal, m 

Superscript referring to maximum applied load 
Superscript referring to minimum applied load 

Range in closure pressure, MPa 

Fiber radius, m 

Stress ratio = S min /S max 

Applied stress range, MPa 
Fiber and matrix volume fractions 

Integration variable along crack from center, mm 
Discrete crack opening displacement range in matrix, Jim 

Continuum crack opening displacement range in composite, |im 

2 

Composite fracture surface energy = / E^, MPa*m 

2 

Matrix fracture surface energy = K m / E m , MPa*m 


x Interfacial frictional shear stress, MPa 

t c Frictional shear stress used to fit the crack growth data, MPa 

t C0( j Frictional shear stress used to fit the crack opening displacement data, MPa 

T^j Frictional shear stress used to fit the debond length data, MPa 

V LT Composite Poisson's ratio 

<Tf, <T m Axial stress in fiber and matrix, MPa 


Introduction 

Fatigue damage progression in advanced titanium matrix composites (TMC) must be properly 
characterized for these materials to be confidently used in a man-rated aircraft. Numerous investigations 
have been conducted on the fatigue damage growth behavior in TMC containing stress concentrations 
(see for example [1-9]). In general, fatigue damage is a complex process which depends on many 
variables including constituent properties, lay-up, fabrication processes, applied loadings, and specimen 
geometries. The dominant mechanisms of fatigue damage in TMC are fiber breakage, matrix cracking, 
and fiber-matrix debonding. Typically, the onset of fiber breakage results in rapid, self-similar damage 
progression and Mode I, catastrophic fracture [1-3]. Under loading conditions where fibers do not 
break, damage progression consists primarily of matrix cracking and fiber-matrix debonding [4-9]. 
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Extensive damage with multiple matrix cracks growing parallel to each other has been observed in TMC 
subjected to maximum loads as low as 17% of the st atic notched strength [5-8], Even though the fibers 
are intact, matrix cracking and fiber-matrix debonding significantly reduced the composite longitudinal 
stiffness and strength [7,9]. In addition, both composite toughness and environmental protection of the 
fibers are lost due to matrix cracking. Thus, matrix cracking is of particular concern. 

Conventional fracture mechanics characterization of Mode I fatigue crack growth behavior is 
accomplished through the relation between the crack growth rate and the stress intensity factor range. 
The most common relation is the power law function proposed by Paris et al. [10]: 

§ = C(AK)" (1) 

where C and n are material constants and AK is the Mode I stress intensity factor range. To apply 
Equation (1) to matrix fatigue cracking in TMC requires the appropriate definitions of AK, C, and n. 
Ideally, the values of C and n should be identical to those for the neat matrix material. Typically, AK is 
equal to the applied stress intensity factor range: 

AK app = AS ViTa (2) 

where AS is the applied far-field stress range. However, matrix cracking in TMC is unlikely to be 
governed solely by AK^^. The stress intensity factor range experienced by the matrix material (AK m ) 

should be a more suitable definition for AK. 

Expressions for AK m for matrix cracks bridged by intact fibers in an unidirectional composite 

were derived in several fiber bridging models [11-15]. As illustrated in Figure 1, these models assume 

that the fiber-matrix interface debonds as matrix cracks progress past the intact fibers. The models also 

assume that a constant, but unknown, frictional shear stress acts over a slip length in the debond region. 

The slip length does not necessarily have to be equal to the debond length. An energy balance approach 
was used by Aveston et al. [1 1] and Budiansky at al. [12] to derive an expression for K m in terms of the 

composite microstructural parameters (constituent moduli, fiber volume fraction, unknown constant 
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frictional shear stress, etc.) under conditions of steady state cracking during monotonic loading. The 
steady state stress intensity factor derived in [1 1,12] is independent of crack length. These two models 
are referred to as steady state fiber bridging (SSFB) models. 

Another class of models, the generalized fiber bridging (GFB) models, combine a continuum 
fracture mechanics analysis and a micromechanics analysis to derive stress intensity factor solutions for 
matrix cracks of arbitrary size. In the GFB models, the constraint due to the intact fibers in the wake of 
the matrix crack is idealized as an unknown closure pressure. The governing equation in terms of the 
unknown closure pressure is obtained by combining the crack opening displacement solutions from the 
continuum fracture mechanics analysis and from the micromechanics analysis. The GFB models 
developed by Marshall et al. [13] and McCartney [ 14 ] were formulated for monotonic loading conditions 
and were modified for fatigue loading conditions by McMeeking and Evans [15]. The GFB models 

differ from each other in the formulations used to relate the continuum fracture mechanics analysis and 
the micromechanics analysis. As the crack length increases, K m asymptotically approaches the steady 

state value given by the SSFB models. 

The objective of this research is to determine the applicability of the GFB models for 
characterizing matrix fatigue crack growth in center notched [0]g SCS-6/TM5-3 and [0]^ SCS-6/Ti-6Al- 

4V laminates. The SCS-6/Ti-15-3 laminates were tested as part of the current study and the SCS-6/Ti- 
6A1-4V laminates were tested by Davidson [16]. MaUix crack initiation and progression were monitored 
and recorded during fatigue loading. The effect of fiber bridging on the matrix stress intensity factor 
range, the crack opening displacement, and the debond length were predicted using the GFB models. 
The unknown frictional shear stress in the GFB models was used as a curve fitting parameter. 
Calculations of the slip length and the crack opening displacement were compared with those measured 
in [16]. In addition, calculations of the slip length were compared to debond lengths measured in this 
study. The Slip lengths in the GFB models were assumed to be equivalent to the debond lengths in this 
study. 
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Materials and Test Procedures 

Materials and Specimens 

The material tested in this study, designated SCS-6/TM5-3, is a titanium matrix composite 

reinforced with continuous silicon-carbide fibers. The composition of the titanium alloy is Ti-15V-3Cr- 

3Al-3Sn. The composite laminates were fabricated by Textron by hot isostatic pressing (HIPing) Ti- 1 c 

3 foils between unidirectional tapes of silicon-carbide (SCS-6) libers having a diameter of 0.14 mm. 
The [0]g laminates had a fiber volume fraction Vj- of 0.33 and were in the as-fabricated condition. Table 

1 lists the material properties for the SCS-6/Ti-15-3 laminates. 

Two SCS-6/TM5-3 specimens were cut using a diamond wheel saw into straight- sided coupons 
with the 0° fibers in the loading direction. Each specimen was 152.4-mm long and 1.8-mm thick. Two 
notch length-to-width ratios (2a/W) were used, 0.30 and 0.35. The center notches were made using 
electro-discharge machining (EDM). To make optical observations and replicas, the surface of each 
specimen was polished to obtain a flat and lustrous finish. Aluminum end tabs were bonded on all 
specimens to prevent specimen failure in the grips. 

Test Procedures 

Constant amplitude, tension-tension fatigue tests were conducted under load control with R = 

0.1 (S • /S ) at a frequency of 10 Hz using a closed-loop servo-hydraulic test machine equipped 
min max 

with hydraulic grips. Matrix crack initiation and progression were monitored and recorded in real time 
using a closed-circuit television system (CCTV) having magnification capabilities up to 325X. Testing 
was periodically interrupted when significant increments in crack extension were observed to take 
surface replicas and to examine the specimens surface using a scanning electron microscope (SEM) and 
an optical microscope. In one specimen, the applied stress level was increased after a prescribed number 
of fatigue cycles while in the second specimen, the applied stress was held constant. The loading history 
of both specimens is given in Table 2. After fatigue loading, the surface of Specimen 1 was polished to 
the midplane of the outer ply of fibers. The mechanisms of fiber-matrix interfacial debonding were 
identified and the debond lengths were measured using the SEM. 



Analytical Models 


The GFB models [13-15] incorporate the dominant fatigue failure mechanisms typically observed 
in TMC, namely, matrix cracking and fiber-matrix debonding. Thus, these models appear to be ideally 
suited to predict the matrix crack growth in TMC. The GFB models combine a continuum fracture 
mechanics analysis and a micromechanics analysis to obtain expressions for the crack opening 
displacement and stress intensity factor for matrix cracks bridged by fibers in a unidirectional composite. 
In the continuum fracture mechanics analysis, crack opening displacement is obtained by modeling the 
constraint of the bridging fibers as a closure pressure. The crack opening displacement is related to this 
closure pressure using the micromechanics analysis and discrete-continuum relations. The governing 
equation reduces to a single nonlinear integral equation in terms of the unknown closure pressure. 

In this section, the derivation of the governing equation is described. First, the continuum 
fracture mechanics and micromechanics analyses are presented. Then, the three discrete-continuum 
relations, designated MCE, MC and ME formulated by Marshall et al. [13], McCartney [14], and 
McMeeking and Evans [15], respectively, are reviewed. Finally, these relations are used to combine the 
continuum fracture mechanics analysis and micromechanics analysis to obtain the governing equation. 


Continuum Fracture Mechanics Analysis 

The continuum solution is obtained by superimposing the solutions of a crack subjected to a far- 
field applied stress range, AS, and a crack subjected to a change in closure pressure, AP(x), as shown in 
Figure 1. The origin of the x-y coordinate system is located at the center of the crack. For a composite 
material, the crack opening displacement is reduced to a single integral equation [17]: 


A5 tip = 


4AS 

E 


- V a 2 -x 2 + 


E'Jt 


AP(x)log 


V a 2 -x 2 W a 2 -x 2 

V a 2 -x 2 -V a 2 -x 2 


dx 


(3) 


where a 0 is the initial crack length without fiber bridging, a is the final crack length, and E' for an 
orthotropic material containing a crack normal to the loading direction is [18]: 
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( 4 ) 



The terms in this equation are the composite laminate properties as defined in the Nomenclature. For an 

isotropic material under plane stress conditions, E' is simply equal to the modulus of the material. \i 

several studies [2,4, 13, 14,19, 20|, E' in Equation (3) was replaced by the longitudinal modulus of the 
composite, E^. However, in this study E' is correctly defined by Equation (4). The composite stress 

intensity factor range for the problem in Figure 1 is [ 17]: 



(5) 


The quantities A§ t jp and AK^p are referred to as the continuum crack opening displacement and the 

continuum stress intensity factor range, respectively, for a composite subjected to fatigue loading and 
can be determined by knowing the closure pressure, AP(x). The quantities A5 t jp and AP(x) are related 

to each other using the micromechanics analysis and the discrete-continuum relations described in the 
following sections. 

Micromechanics Analysis 

AA' 

The schematic in Figure 2a illustrates the fiber stress along the crack center line, A a f , as a 

function of the matrix crack opening displacement, A5 m , during fatigue loading. Three regions are 

shown: (1) initial loading, O-A; (2) unloading, A-B; and (3) reloading, B-A. Marshall and Oliver 

[21] derived relations for each of these three regions and made comparisons with results from fiber push- 

out tests of ceramic matrix composites. McMeeking and Evans [15] used a similar force balance 
approach to determine the A5 m - Aa AA relation. The micromechanics analysis provides the vital link 

between the discrete matrix crack opening displacement, A8 |T| , and the fiber stress along the matrix crack 
centerline, Aa^, as shown in the Appendix. The unknown closure pressure and the fiber stress along 
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the matrix crack centerline, Ac^, are related through the fiber volume traction [13-15]: 


AP = v f (AcJf A ) 


( 6 ) 


This equation is valid if the fibers are closely spaced relative to the crack length. Using Equation (6), the 
matrix crack opening displacement is related to the unknown closure pressure: 


AS m = XAP 2 


(7) 


where: 


X = 


r 

4tv^EjT1 


( 8 ) 


and: 



(9) 


Equations (6) and (7) are valid in the wake regions away from the crack-tip; however, the singularity 
fields in the crack-tip vicinity are not included in the micromechanics analysis. Consequently, Equations 
(6) and (7) erroneously suggest that the fiber stress vanishes at the crack-tip. 


Discrete-Continuum Relations 

The continuum fracture mechanics analysis and the micromechanics analysis are combined using 
the following two discrete-continuum relations: (1) A5 m - A5 lip , and; (2) AKjp - AK^. Differences 

among the GFB models [12-14] arise in defining these two discrete-continuum relations as discussed 

next 
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Crack Opening Displacements 

Marshall ct al. [13] and McMeeking and Evans [ 1 5] assumed that the discrete and the continuum 
crack opening displacements were equal: 

**■, - A*up (10) 

McCartney [14] reported that A5^p is actually the change in displacement within the slip region (change 
in displacement between lines AA' and BB\ Figures 2b and 2c) and should be related to A8 m as 

follows: 

(,1) 

where r\ is defined in Equation (9). Both relations are examined in the present study. 


Stress Intensity Factor 

In order to use a criterion for matrix cracking, a relation between the matrix and continuum stress 
intensity factors must be established. Marshall ct al. [13] related AK m to AK^p using the ratio of the 

matrix modulus to composite modulus: 




( 12 ) 


This equation assumes that the near-tip strains in the composite and in the matrix are compatible. In a 
different approach, McCartney [14] used an energy balance to relate AK m to AK tip by assuming: 

Y = v Y (13) 

*c m'm 

2 2 
where y c = K tip^L * s the ^ racturc su, f ace energy of the composite and y m = is the fracture 
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surface energy of the matrix. Rewriting Equation (13) in terms of stress intensity factor ranges: 

4K - - V5T^ a4) 

Finally, it was suggested by McMeeking and Evans [15]: 

AK m - AK t,p (15) 

All three stress intensity factor relations, Equations (12), (14), and (15), are examined in this study. 

In the subsequent sections, the discrete-continuum relations derived by Marshall et al. [13] (A5 m 

E 

= A5 Up andAK m = ^ AK Up ), McCartney [14] (A8 m = A5 tip Ti and AK m = 

L 

and McMeeking and Evans [15] (A5 m = A8 t j p and AK m = AK^p) are referred to as the MCE, MC, 
and ME relations, respectively. 


V 


m 


c AK tJn ) 
v m E L l, P 


Governing Equation 

The governing crack opening displacement equation is obtained by substituting either Equation 
(10) or (1 1) into (3). After normalizing, the governing equation becomes: 


P(u) 


- mVTV = - ^ 


P(v)log 


V 1-U 2 W 1-v 2 




dv 


(16) 


where: P = AP/AS 


u = x/a 
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For either the MCE or ME displacement relation (Equation (10)), the term p is: 

4a 

“"SS (17) 

and for the MC displacement relation (Equation (11)), the term p is: 

11 = iSP < 18 > 

Equation (16) is a nonlinear integral equation that is solved numerically using an iterative procedure 
similar to that outlined in [14]. Once the normalized pressure is known, the matrix crack opening 
displacements are found from either the MCE and ME displacement relations (Equation (10)) or from the 
MC relation (Equation (11)). Then the matrix stress intensity factor ranges are obtained from the MCE, 
MC, and ME stress intensity factor relations. Equations (12), (14) and (15), respectively. 

Results and Discussions 

In this section, the frictional shear stress in the GFB model was used as a curve fitting parameter 
to matrix crack growth data, debond length data, and crack opening displacement data. The three 
discrete-continuum relations (i.e. the MCE, MC and ME relations) were used in the analysis. For each 
TMC system, the experimental results are presented first, followed by the curve fit results using the 
GFB model. Interpretations of the lrictional shear stress term in the GFB model are discussed as well as 
the differences in the frictional shear stress values used to fit the data. Finally, the limitations and 
advantages of the GFB model are discussed. 

SCS-6/Ti-15-3 
Matrix Crack Growth Data 

The average cumulative crack extension as a function of the number of cycles is shown in Figure 
3 for both SCS-6/Ti- 15-3 specimens tested. In Specimen 1 the applied stress level was increased after a 
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prescribed number of fatigue cycles and in Specimen 2 the applied stress was held constant, as shown in 

Table 2. In Specimen 1, the crack growth rate increased as the applied stress level increased. For both 

specimens, the crack growth rate decreased as the crack length increased, as shown in Figure 4. 

A typical matrix crack for Specimen 1 subjected to S max = 300 MPa is shown in Figure 5. The 

crack shown in the first photograph was developed during the previous load history. As shown in this 

figure, the matrix crack progressed from the notch-tip in a macroscopically self-similar manner. It was 

assumed t hat the fibers in the wake of the matrix crack were intact since no jumps in crack opening 

displacement (COD) were observed and no audible levels of acoustic emission were heard (fiber 

breakage is associated with sudden increments in COD and high amplitude acoustic emission events [1]). 
Naik and Johnson [6] also did not observe any fiber breaks in [0]g SCS-6/Ti-15-3 specimens containing 

double edge notches subjected to similar loading conditions when the outer layer of matrix material was 
etched away. The intact fibers in the wake of the matrix crack effectively reduced the crack growth rate 
as the crack length increased. 

The matrix fatigue crack growth behavior in the composite was first characterized using the neat 
matrix material properties and the applied stress intensity factor range (AK ), Equation (2) calculated 

a rr 

with standard data reduction procedures (ASTM Standard E647). The crack growth rate (da/dN) as a 
function of AK,^ f° r the Ti-15-3 sheet material [22] and the two SCS-6/Ti-15-3 specimens is shown in 

Figure 6. As shown in this figure, using AK a pp in the composite significantly overestimates the actual 

stress intensity factor range governing matrix crack growth. In addition, there is a trend towards a 
negative slope in the da/dN - AK a pp curves for the two composite specimens. For a given AK a pp, the 

da/dN is reduced by approximately one order of magnitude due mainly to the bridging fibers compared to 
the neat matrix results. 

The GFB models were then used to determine the matrix stress intensity factor range (AK m ) in 

the two specimens tested. Using the unknown frictional shear stress as a curve-fitting parameter, the 
matrix crack growth data was collapsed onto that of the Ti-15-3 sheet material as shown in Figures 7 
through 9. The crack growth data was reduced in Figures 7, 8 and 9 using the MCE, MC and ME 

relations, respectively. Each figure indicates the discrete-continuum relation used as well as the values of 

% the frictional shear stress required to fit the composite crack growth data. As shown in these 
eg 
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figures, the calculated stress intensity factor was reduced when bridging fibers were modeled. In 
addition, x C g typically increased as the applied stress increased. There was a wide range in the values of 

x C g needed to fit the experimental data, as shown in these figures. 

For the specimen subjected to increasing stress levels (Specimen 1), an excellent fit was obtained 
for each stress level using the value of x C g shown in Figures 7 through 9. However, for the specimen 

subjected to a constant applied load (Specimen 2), the fit was not as good for the MC and ME relations, 
Figures 8 and 9 respectively. In Specimen 1, the crack extensions that occurred during each load level 
were short compared to that in Specimen 2 (see Table 2). For Specimen 2, a better fit would have been 
obtained if larger values of x were used at the higher crack growth rates (corresponding to small crack 
lengths). This would imply that the frictional shear stress decreases as the crack length increases. 

Debond Length Data 

In order to study the mechanisms of fiber-matrix debonding, the surface of Specimen 1 near the 
notch-tip was polished to the midplane of the outer ply and was examined under an SEM. During 
consolidation, a brittle reaction zone developed between the matrix and the outer of the two carbon rich 
coatings of the SiC fiber as shown in Figure 10. In the interface of the first intact fiber, the reaction zone 
near the matrix crack had a rubble-like appearance indicating fracture of the brittle reaction zone. Further 
away from the matrix crack along the fiber (approximately 0.5 mm away), the reaction zone had a more 
uniform, intact appearance. However, cracking was now observed between the carbon rich layers. 
Cracking in the carbon rich layers extended approximately 1.5 mm along the first intact fiber. The 
transition from fracture of the brittle reaction zone to cracking between the carbon rich coatings occurred 
approximately 0.42 mm from the crack centerline as shown in Figure 11. In the other bridging fibers, 
cracking of the carbon rich coating was less prevalent and the debond length was composed mainly of 
fracture of the brittle reaction zone. 

Measurements of the debond length for each fiber in the bridged region are shown for Specimen 
1 in Figure 12. The ordinate axis is the distance measured from the first intact fiber. The entire interface 
damage process zone (failure of reaction zone and cracking in the carbon rich coatings) was assumed to 
be the debond length. The precise end point of the debond region was difficult to identify; thus, the 
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results presented in Figure 12 are approximate. As shown in Figure 12, the debond lengths for the 
fibers in the bridged region decreased as the distance from the first intact fiber increased. 

In the GFB models, the slip length is defined as the region over which the frictional shear stress 
acts. In the present study, the slip length is assumed to be equal to the debond length, and the measurer 4 

debond lengths are compared to calculated slip lengths. For Specimen 1 the calculated slip lengths are 
presented in Figures 13 and 14 assuming A8 m = A5 t jp (MCE and ME relations) and A5 m = A5^pT] 

(MC relation), respectively. Here, calculations of the slip length are made using several values of 
frictional shear stress including the values used to fit the crack growth data, x C g, and the debond length 

data, Xjj. In Figure 13, the two values of x C g corresponding to the MCE and ME stress intensity factor 

relations are labeled. In general, as x increased, the slip length decreased, where the slip length is 
proportional to The best agreement between the calculations and experiments was obtained for x^j 

= 40 MPa. For all discrete-continuum relations, the debond lengths were overestimated using the 
corresponding values of x C g as shown in Figures 13 and 14. The best agreement between x^j and x C g 

was obtained assuming the ME relations [15] (x C g = 20.43 MPa and x^j = 40.00 MPa, Figure 13). 

As shown in Figures 1 3 and 1 4, a better agreement between the measured debond lengths and the 
calculated slip lengths for the first two fibers in the bridging region would be obtained if smaller values 
of x were used. In fact as was shown in Figure 10, the debond length for the first intact fiber consisted 
mainly of cracking in the carbon rich layers. Although not shown, in the other bridging fibers 
debonding consisted mainly of fracture of the brittle reaction zone. The crack surfaces in the carbon 
layers were smoother than the fracture surface in the brittle reaction zone. Hence, in the first intact fiber, 
it is possible that the overall shear stress acting on the debond surfaces would be lower . 

SCS-6/Ti-6Al-4V 
Matrix Crack Growth Data 

The effects of the frictional shear suess on the matrix crack growth data were determined for the 
[0] 4 SCS-6/Ti-6Al-4V laminate tested in [16]. The reduced matrix crack growth data is shown in 

Figures 15 , 16, and 17 using the MCE, MC and ME relations, respectively. Here, predictions of the 
matrix stress intensity factor range are made using several values of frictional shear stress including 
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values used to Fit the crack growth data , x C g, the debond length data, x^j, and the crack opening 
displacement data, x C0( j. As shown in these figures, the applied stress intensity factor range 

overestimates that of the matrix and, therefore, cannot be used to characterize the matrix crack growth 

behavior. Using the GFB models, the crack growth data for the in situ matrix material and the neat 
matrix material were forced to coincide by varying x . The best agreements between experiment and 

predictions were obtained for x C g = 0.9 MPa, 1.0 MPa, and 4.5 MPa for Figures 15, 16 and 17, 

respectively. Calculations of the matrix stress intensity factor range was sensitive to the value of x used. 

As x increased, the matrix stress intensity factor decreased, where the matrix stress intensity factor range 

was proportional to 1 /^fx. An increase in X results in an increase in the closure pressure. Consequently. 

both the crack opening displacement and the stress intensity factor were reduced when x increased. 

Using the MCE and MC relations, the matrix stress intensity factor range was underestimated 
using x^j and t C0( j as shown in Figures 15 and 16, respectively. The best agreement among the 

calculations using x^ , x CQ( j and x £ g was obtained using the ME relations as shown in Figure 17. In this 

figure, the experimental values of the crack growth data are bounded by the predictions made using x^j 

(lower bound) and x co< j (upper bound). The value of AK m varied approximately ±3 MPaVm between 

these bounds. 

Debond Length Data 

The slip lengths shown in Figures 18 and 19 were calculated assuming A5 m = A5^p (MCE and 
ME relations) and A5 m = AS^ri (MC relation), respectively. In these figures, the ordinate axis is the 

length measured from the first intact fiber. Calculations of the slip length are made using several values 
of frictional shear stress including x C g, x^, and x coc j. Comparisons are made with the measured debond 

lengths. In general, the calculated slip length decreased as x increased, where the calculated slip length 

was proportional to 1/V^ In Figures 18 and 19, the best fit between the experiments and calculations 
was obtained for x dl = 12.5 MPa. In Figure 18, a reasonably good fit was obtained using a single 

value of the frictional shear stress. However, in Figure 19, a better fit would be obtained if a larger 
value of Xjjwere used for the first two intact fibers in the bridging region. Thus, x^ is apparently a 

function of the distance from the first intact fiber for this material as well as the SCS-6/TI-15-3 tested. 
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For all discrete-continuum relations, the debond lengths were overestimated using the 
corresponding values of x cg , as shown in Figures 18 and 19. In Figure 18, two values of x QO are 

labeled corresponding to the MCE and ME stress intensity factor range relations. The best agreement 
between x d j and x cg was found assuming the ME relations (x £g = 4.5 MPa and x d j = 12.5 MPa, Figure 

18). In this case, the calculated debond lengths using these values of x cg were approximately twice as 

large as the experimental measurements. Figure 18. For all three discrete-continuum relations, the 
debond length was severely overestimated using x cod as shown in Figures 18 and 19. 


Crack Opening Displacement Data 

The crack opening displacement shown in Figures 20 and 21 were calculated assuming AS m = 
^tip (MCE and ME relations) and A5 m = A5 d pT| (MC relation), respectively. In these Figures, the 

ordinate axis is the length measured from the First intact fiber. Calculations of the crack opening 
displacement are made using several values of frictional shear stress including x cg , x d j, and x cod . In 

Figure 20, two values of x £g are labeled corresponding to the MCE and ME stress intensity factor 

relations. In general, the calculated crack opening displacements decreased as X increased. The 

calculated crack opening displacement was proportional to 1/x. The best Fit between experiments and 
calculations was found for x C()d = 1.5 MPa and x cod = 3.0 MPa. as shown in Figure 20 and 21, 

respectively. In both Figures, a better Fit would have been obtained if larger values of x cod were used for 

the first four intact Fibers in the bridging region. For all discrete-continuum relations, the crack opening 
displacement data was underestimated using x d j. For the MCE and MC relations, the crack opening 

displacement data was bounded by the calculations made using x d j and x cg . Figures 20 and 21, 
respectively. 

Interpretation of Frictional Shear Stress 

The unknown constant frictional shear stress used in the fiber bridging models is a critical 
parameter. As shown previously in Figure 10, cracking of the reaction zone occurs near the matrix 
crack, but further along the fiber-matrix interface, a transition in the debonding mechanism was observed 
and cracking between the carbon layers was found. The value of x should be considerably different in 
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ihese two debonded regions. The value ot X used in the GFB models can only represent an average of 

the actual shear stress distribution acting on the debonded surfaces. 

The actual frictional shear stress is a difficult cjuanlily to measure and, thus, it has been frequently 
used as a curve fitting parameter. However, x is related to both the debond length and the crack opening 
displacements, which can be directly measured. The inlerfacial frictional shear stress can only have some 
physical significance if the same value of x can be used to accurately predict the matrix crack driving 

force, the debond length, and the crack opening displacement. The values of x used to fit the crack 
growth data (x cg ), the debond length data (x dJ ), and the crack opening displacement data (x cod ) are listed 

in Table 3 according to the discrete-continuum relation used in the calculation. No crack opening 

displacement data was measured for the SCS-6/Ti- 15-3 laminates tested here. As indicated in this table, 
there is a large difference in the values of x cg , x cod and x d , even within a specific discrete-continuum 

relation. Thus, quantifying the frictional shear stress using the GFB models does not appear to be a 
valid approach. 

The ME relations, A5 m = A5 lip and AK m = AK tip , produced the best fit for the data sets in 

this study and in [16] with the least variation in the frictional shear stress. Using the ME relations, all 

values of X were between 4 and 360 MPa, the range of x reported in [2,15,20,23,24,25,26] for both 
SCS-6/Ti-15-3 and SCS-6/Ti-6Al-4V, except for the value of x cod (1.5 MPa) for the SCS-6/Ti-6Al-4V 

specimen. From fiber push-out tests, the reported range of values of x was from 30 to 360 MPa 
[20,25,26] which is an upper limit to the values calculated in this study. The large variation among the 
reported values of frictional shear stress measured from fiber push-out tests is probably due to 
differences in test methodologies, composite fabrication procedures, etc. In addition, the measured 
frictional shear stress values from the fiber push-out tests may not accurately represent the actual shear 
stress acting along the debond region in a bridging fiber. Thus, comparisons with values of x calculated 

in this study may not be appropriate. 

Advantages and Limitations of Fiber Bridg ing Models 

The evaluated fiber bridging models do capture the most essential features of the fatigue damage 
progression in the materials studied, i.e., both matrix cracking and fiber-matrix debonding are modeled. 
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Thus, an expression for the stress intensity factor range driving the observed matrix cracks was 
obtained. In addition, these models are relatively simple to program and can be used to determine the 
severity of damage based on the fiber stress or the matrix stress intensity factor. The matrix fatigue 
cracking can be effectively characterized without using complicated numerical methods, such as finite 
element analyses, which can be quite cumbersome. 

The GFB models were useful for the parametric studies conducted here and for interpreting 
experimental results. In another useful application, the effects of fiber-matrix interface strength on the 
fatigue behavior can be studied by varying x. The el feci ol the frictional shear stresses on the stress 

intensity factor range and the stress in the first intact fiber is shown in Figure 22. As x increases 
(increasing fiber-matrix bond strength) the fiber stress increases while AK m decreases. Thus, the fibers 

would tend to fracture in a composite with a strong fiber-matrix interface, whereas, matrix cracks would 
more likely propagate in a composite with a weak liber-matrix interface. This trend agrees with 
observations made by Naik et al. [281 on the effects of the fiber-matrix interface strength on the fracture 
and fatigue properties of [0/90/0] SCS-6/Ti-15-3 laminates. 

There are several limitations to the fiber bridging models worth noting. The fibers bridging the 
matrix cracks are idealized as a continuous closure pressure. The micromechanical analysis relating the 
closure pressure to crack opening displacement is based on a simplified one-dimensional analysis of a 
fiber in the wake of the matrix crack and does not consider the crack-tip mechanics. The complex 
micromechanical details at the crack tip are not modeled. Consequently, the fiber stresses calculated 
using the fiber bridging models vanish at the crack-tip, which is unrealistic. In addition, the fiber 
bridging models cannot take into account the three dimensional effect of the crack front bowing around 
the fibers as discussed and modeled by Bower and Ortiz [27]. This toughening mechanism, termed 
"crack trapping", can considerably reduce the crack driving force. It was shown in [27] that crack 
trapping can double the toughness of a brittle matrix composite without fiber bridging. Finally, 
modeling the frictional shear stress as a constant is approximate at best. As discussed earlier, the 
variation in debonding mechanisms along the interface would certainly yield differences in X. In 
addition, the large variation in the values of x used in fitting the various data sets is quit disturbing. The 
frictional shear stress is apparently a function of crack length, applied stress, distance from the first intact 
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liber, and distance along the debond length, and, thus, is not a material constant. 

The predictive capabilities of the GFB models are very suspect due to the dependency of x on so 
many factors. However, the GFB models may provide a frame work for a crack growth prediction 
methodology. This may be accomplished by incorporating the crack-tip singularity fields in the 
micromechanics analysis and reevaluating the discrete-continuum relations which appear to be over- 
simplified. In addition, a value of the frictional shear stress should be established through experiments 
and used in the analysis since it resembles a material parameter. A new fitting parameter can then be 
introduced in the model to include mechanisms not currently being considered such as crack trapping. 

Concluding Remarks 

In this study, the application of the generalized fiber bridging (GFB) models to analyze matrix 
crack growth in [0]g SCS-6/Ti-15-3 and [0]^ SCS-6/Ti-6Al-4V laminates was investigated. 

Experimental observations revealed that the fatigue damage in [0]g SCS-6/TM5-3 laminates consisted 

primarily of Mode I matrix cracks that initiated and grew from notch-tips. Fibers were intact in the 
wake of the matrix cracks. In addition, fiber-matrix debonding in the form of fracture of the brittle 
reaction zone and cracking between the two carbon rich fiber coatings was found. Under a constant 
applied stress range, crack growth rate decreased as the matrix crack length increased. The fibers in the 
wake of the matrix cracks effectively reduced the stress intensity factor. The applied stress intensity 
factor range did not characterize the fatigue crack growth of matrix cracks which were bridged by fibers. 

The effect of the fibers bridging the matrix crack was studied using the GFB models. In these 
models, fiber-matrix debonding is modeled as the matrix crack progresses past the fibers. An unknown 
constant frictional shear stress acts over the debond or slip length. These models derive the matrix stress 
intensity factor by combining a continuum fracture mechanics analysis and a micromechanics analysis 

using discrete-continuum relations. In the continuum fracture mechanics analysis, expressions for the 
composite stress intensity factor range (AK tip ) and the composite crack opening displacement (A8 tip ) are 

derived. The micromechanics analysis provides the solution to the matrix crack opening displacement 
(A8 m ) in terms of the unknown closure pressure and the composite micromechanical parameters. 

Discrete-continuum relations are used to relate AK m to AK t j p and A8 m to A8 t - p . 


19 



The effect of fiber bridging on the matrix stress intensity factor range, the crack opening 
displacement, and the debond length were studied using the GFB models. The GFB models were 
exercised using the unknown constant frictional shear stress, x, as a curve fitting parameter to available 

data. In general, as x increased, the calculated matrix stress intensity factor, crack opening displacement, 
and slip length all decreased. By assuming AK m = AK tip and A8 m = A5 tip , the GFB models yielded 

the most accurate correlations based on the smallest scatter in Irictional shear stress. However, there was 
a large difference in the values of x used to correlate the data even within a specified discrete-continuum 
relation. The value of x depended on the crack length, applied stress level, and distance from the first 
intact fiber and is not a material property. 

There are several shortcomings of the GFB models. The fiber bridging models cannot account 
for the three-dimensional effects such as crack front bowing, which can considerably reduce the matrix 
stress intensity factor. The GFB models assume the frictional shear stress acting along the debonded 
fiber-matrix interface is a constant. The different debonding mechanisms observed along the interface 
should, however, result in a variation in the frictional shear stress. Finally, the predictive capabilities of 
the GFB models are questionable due to the dependency of x on many factors. 

In spite these limitations, the GFB models do incorporate the major modes of fatigue damage 
(matrix cracking and fiber-matrix debonding) and yield expressions for the stress intensity factors that 
characterize the matrix crack driving force. The GFB models provide an efficient and relatively simple 
engineering approximation to conduct parametric analysis using the composite micromechanical 
variables. 
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Appendix: Micromechanics of Fiber Bridging 

The relation between the discrete crack opening displacement and the fiber stress along the crack 
centerline is derived in detail for each region of fatigue loading shown in Figure 2a. A free body diagram 
of a fiber in the wake of the matrix crack for each loading phase is shown in Figures 2b and 2c. 


Loading: Q-A 

During the initial loading along O-A (Figure 2a), the frictional shear stress (x) is constant in the 
slip region (i ) and opposes the fiber stress, Figure 2b. This free body diagram is constructed assuming 
the maximum load is just reached (Point A in Figure 2a). Along line BB' (end of slip region), the strains 
in the fiber and the matrix are equal: 


<°r" 


max 


( BB^max 

v m ' 


"m 


(Al) 


where (of B ') max and (o“ B ' ) m ax arc ^ fiber and matrix stresses along line BB' at the maximum applied 
load (Point A in Figure 2a). Overall equilibrium in the matrix in the slip region requires: 

(°m B ) maXA m * 2kxv!1 = 0 (A2) 


where A m is the cross-sectional area of the matrix. Equilibrium in the fiber in the slip region requires: 
(o“' + 2lttr . i = «jf A ') maX Af (A3) 

where (cr^ ) max is the fiber stress along the matrix crack centerline AA' and A m is the cross-sectional 

area of the fiber. Combining Equations (Al), (A2), and (A3), the fiber stress along the crack centerline 
AA' at the maximum load is: 
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The total extension of the matrix in the slip region is: 
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Combining Equations (Al), (A2), (A6), and (A7), an expression for the slip length, SL, is obtained in 
terms of the crack opening displacement at the maximum load, 5 max : 


SL = 


.lC rE f 

V 2xn 


(A8) 


Substituting Equation (A8) into (A4) yields the relation between the fiber stress along the crack centerline 
and the discrete crack opening displacement at the maximum applied load: 
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Unloading! A-B 

Upon unloading, the frictional shear stress reverses direction within a length, i R , in the slip 

region. The free body diagram. Figure 2c, is constructed assuming the minimum load is just reached 
(Point B in Figure 2a). The procedure to determine the relation between the fiber stress and discrete 
crack opening displacement is similar to that at maximum load. Along line BB' (end of slip region), the 
strains in the fiber and the matrix are equal: 



(A10) 


Overall equilibrium in the matrix in the slip region requires: 
(C) mmA m ’ 27txr(i-i R )+ 2fttrl R s 0 


(All) 


and in the fiber: 

(aj? B ) m * n Af + 27tir(i - -6 r ) - 2n:xri R = (a AA ) min Aj- (A12) 

Combining Equations (A 10), (All), and (A 12), the fiber stress along the crack centerline at the 
minimum load is: 
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(A 13) 


Substituting Equation (A4) into (A13): 
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(A 14) 




where Acj AA = (a AA ) max - (a AA ) min . The total extension of the matrix in the slip region at the 
minimum load is: 
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(A 15) 


which reduces to: 
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The extension of the fiber at the minimum load is: 
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which reduces to: 
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(A 18) 
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Combining Equations (A 10), (All), (A 16), and (A 18): 


min _ 2rp jo 2 a 2, 
°m - re J* ) 
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(A19) 


By substituting Equation (A8) into (A 19), the length of the shear stress reversal at the minimum load is: 


^ r= VS 


A8 


m 


(A20) 


max .min 

where A8 m = 8 m - 5 m . Combining Equations (A 14) and (A20) yields the relation between the 

change in fiber stress along the crack centerline and the change in discrete crack opening displacement 
from points A to B in Figure 2a: 


AS = _J_ {Aaf A ] 2 
m 4tE f nl f J 


(A21) 


Reloading: B-A 

When reloading back to A, the constant frictional shear stress in <0^ reverses direction. Upon 

reaching A, the shear stress within the entire slip length, J0, is in the same direction. The free body 
diagram in this case is the same as that shown for the original loading case, Figure 2b. Consequently, 
the discrete crack opening displacement when reloaded to Point A is identical to Equation (A9). 
Moreover, the change in displacement from Point B to A in Figure 2a is the same as that from Point A to 
B (Equation (A21)). 
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fable 1. Material Properties: 



Property 

SCS-6/TM5-3 

SCS-6/Ti-6Al-4V [16] 

Matrix: 

E m (GPa) 

92.4 

110.0 


V 

0.34 

0.34 


v m 




C (Vm/MPa*cycle) 

6.49E-1 1 

4.00E-1 1 


n 

2.72 

3.00 

Fiber: 

E f (GPa) 

400.0 

400.0 


v f 

0.25 

0.25 

Composite: 

Laminate 

[0] 8 

[0] 4 


v f 

0.33 

0.42 


E l (GPa) a 

192.2 

231.8 


Ey (GPa) a 

130.1 

165.5 


G lt (GPa) a 

45.9 

59.6 


v a 

V LT 

0.32 

0.30 


a Predicted based on constituent properties [29] 
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Table 2. Specimen Loading History: 


Specimen 

a Q (Initial) 

Number 

(mm) 

1 

3.048 

1 

3.254 

1 

3.572 

1 

3.745 

2 

3.316 

2 

3.375 


aj- (Final) 

Applied Stress 

(mm) 

W MPa > 

3.254 

65 

3.572 

120 

3.745 

200 

4.281 

300 

3.375 

65 

5.081 

325 


Number 
of Cycles 


250,000 (precracking) 

150.00 

40.000 

50.000 


250,000 (precracking) 
250,000 
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30 


1.50 MPa 3.00 MPa x .=1.50 MPa 
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Figure 1. Continuum fracture mechanics approximation of matrix crack growth process during fatigue loading in TMC. Fibers 
bridging crack are idealized as a closure pressure. Governing equation is formulated by superposition of a crack 
subjected to a far-field applied stress and a crack subjected to a closure pressure. 


Fiber Stress 
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Figure 2. Micromechanics of fiber bridging: (a) Fiber stress as a function of crack opening displacement; 
(b) Free body diagram of fiber in wake of crack at maximum applied load; and (c) Free body 
diagram of fiber in wake of crack at minimum applied load. 


SCS-6/TI-15-3 
Lay-Up: [0^ 

v, = 0.33 
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Figure 3. Average cumulative crack extension as a function of the number of fatigue cycles. As the number of cycles 
increases, the rate of crack growth decreases due to fibers bridging the matrix crack. 


SCS-6/TM5-3 
Lay-Up: [0],, 
v# = 0.33 
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Figure 5. Development of notch-tip matrix cracks during fatigue loading, = 0.33. 
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Figure 6. Crack growth rate as a function of the applied stress intensity factor range. The applied stress intensity factor 
range overestimates the stress intensity factor range in the matrix and, therefore, cannot be used to characterize 
matrix crack growth. 
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Figure 7. Crack growth rate as a function of stress intensity factor range for given values of applied stress. The stress 
intensity factor range in the matrix was calculated using the MCE discrete-continuum relations for the given 

values of T . 
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Figure 8. Crack growth rate as a function of stress intensity factor range for given values of applied stress. The stress 
intensity factor range in the matrix was calculated using the MC discrete-continuum relations for the given 

values of x. 
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Figure 9. Crack growth rate as a function of stress intensity factor range for given values of applied stress. The stress 
intensity factor range in the matrix was calculated using the ME discrete-continuum relations for the given 

values of % . 




CO 

a> 

o 

>» 

O 

o 

o 

^ N 

— ox 
n _ 

O CM O 


it ii 

X Z 


II 


00 


O 

«— i (Q 

_ ji > 

CO «e 

in 5 

V « § 
p « § 
co n 

cn » 8 

O ^ E 

<0 :T(/) 



40 


cracking between the carbon rich coatings. 



41 


Figure 11. Mechanisms of fiber-matrix debonding along the first intact fiber. Transistion of fiber-matrix debonding 
mechanisms occurred at approximately 0.42 mm from the matrix crack. 
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Figure 15. Crack growth rate as a function of stress intensity factor range. The applied stress intensity factor range overestimates 
the AK governing matrix crack growth. Using the MCE discrete-continuum relations, AK m was calculated for given 
values of x. Best agreement between experiments and calculations was obtained for x C g- 0.9 MPa. 
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AK (MPa VnT ) 

Figure 16. Crack growth rate as a function of stress intensity factor range. The applied stress intensity factor range overestimates 
the AK governing matrix crack growth. Using the MC discrete-continuum relations, AK m was calculated for given 
values of t. Best agreement between experiments and calculations was obtained for t ra = 1.0 MPa. 
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Figure 17. Crack growth rate as a function of stress intensity factor range. The applied stress intensity factor range overestimates 
the AK governing matrix crack growth. Using the ME discrete-continuum relations, AKj^was calculated lor given 

values of x. Best agreement between experiments and calculations was obtained for x C g= 4.5 MPa. 
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Figure 18. Slip length calculated using the MCE and ME displacement relations for given values of t. As x increases, the slip 
length decreases. Best agreement between experiments and calculations was obtained for x H1 = 12.5 MPa. 
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Figure 20. Crack opening displacement calculated using the MCE and ME displacement relation for given values of x. As x 
increases, the crack opening displacement decreases. Best agreement between experiments and calculations was 
obtained for ^ C0( j= 1.5 MPa. 
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Distance From First Intact Fiber (mm) 

Figure 21. Crack opening displacement calculated using the MC discrete-continuum relations for given values of x. As x 
increases, the crack opening displacement decreases. Best agreement between experiments and calculations v. 
obtained for x„ nA = 3.0 MPa. 
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Figure 22. Stress intensity factor range in matrix and stress in first intact fiber as a function of the interfacial frictional shear. As 

x increases, the stress in the first intact fiber increases, thus, promoting fiber breakage. As x decreases, the matrix 
stress intensity factor range increases, thus, promoting matrix cracking. 
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